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ABSTRACT 

Recent numerical simulations of coalescences of highly spinning massive black hole 
binaries (MBHBs) suggest that the remnant can suffer a recoil velocity of the order 
of few thousands km/s. We study here, by means of dedicated simulations of black 
holes build-up, how such extreme recoils could affect the cosmological coalescence 
rate of MBHBs, placing a robust lower limit for the predicted number of gravitational 
wave (GW) sources detectable by future space-borne missions (such as LISA). We 
consider two main routes for black hole formation: one where seeds are light remnants 
of Population III stars (~ 10^ M©), and one where seeds are much heavier (>10^ Mq), 
formed via the direct gas collapse in primordial nuclear disks. We find that extreme 
recoil velocities do not compromise the efficient MBHB detection by LISA. If seeds are 
already massive and/or relatively rare, the detection rate is reduced by only ~ 15%. 
The number of detections drops substantially (by ^ 60%) if seeds are instead light 
and abundant, but in this case the number of predicted coalescences is so high that 
at least ^10 sources in a three year observation are guaranteed. 
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1 INTRODUCTION 

• Massive black hole (MBH) binaries (MBHBs) are among the 
' primary candidate sources of gravitational waves (GWs) at 
, mHz frequencies, the range probed by the space-based Laser 

• Interferometer Space Antenna {LISA, Bender et al. 1998). 
Today, MBHs are ubiquitous in the nuclei of nearby galax- 
ies (see, e.g., Magorrian et al. 1998). If MBHs were also 
common in the past, and if their host galaxies experienced 
multiple mergers during their lifetime, as dictated by popu- 
lar cold dark matter hierarchical cosmologies, then MBHBs 
inevitably formed in large numbers during cosmic history. 
Provided MBHBs do not "stall", their inspiral driven by 
radiation reaction follows the merger of galaxies and proto- 
galactic structures at high redshifts. MBHBs coalescing in 
less than a Hubble time would give origin to the loudest 
GW signals in the Universe, and a low-frequency detector 
like LISA will be sensitive to GWs from binaries with total 
masses in the range 10'^ — 10^ M© out to 2: ~ 20 (Hughes 
2002). 

The formation and evolution of MBHs has been in- 
vestigated recently by several groups in the framework of 
hierarchical clustering cosmology (e.g. Menou, Haiman & 
Narayanan 2001, Volonteri Haardt & Madau 2003, Koushi- 
appas. Bullock & Dekel 2004). The inferred LISA detection 
rate, ranging from a few to a few hundred per year, were de- 
rived in a number of papers (Jaffe & Backer 2003, Wyithe & 



Loeb 2003, Sesana et al. 2004, Sesana et al. 2005, Enoki et al. 
2004, Rhook & Wyithe 2005). More recently Sesana Volon- 
teri & Haardt (2007) investigated the imprint of massive 
black hole formation models on the expected MBHB coales- 
cence rate, finding that at least ~ 10 (considering a model 
that marginally reproduces the observational constrains and 
that can be taken as a robust lower bound) sources should 
be safely regarded as observable by LISA, assuming a 3 year 
lifetime mission. 

GWs emitted during the final plunge of the binary, carry 
away a net linear momentum, causing a recoil of the MBHB 
center of mass in the opposite direction (Redmount & Rees 
1989). This GW recoil could have interesting astrophysi- 
cal effects, since many coalescence remnants can be ejected 
from their host galaxies and dark matter halos (e.g. Madau 
et al. 2004, Merritt et al. 2004, Micic Abel & Sigurdsson 
2006). This justifies the increasing effort to obtain accurate 
estimates of the recoil velocity. In the case of non spinning 
black holes, the latest analytical (e.g. Favata Hughes & Holz 
2004, Blanchet Qusailah & Will 2005, Damour & Gopaku- 
mar 2006) and numerical (Baker et al. 2006, Gonzalez et 
al. 2007) approaches are now both converging to maximum 
recoil velocities Vr in the range 100-250 km/s for binaries 
with mass ratio q = m2/m\ ~ 0.4 (m2 < mi are the masses 
of the two binary members). The expected values are only 
slightly higher if the binary is eccentric; Sopuerta Yunes and 
Laguna (2007) found oc (1 + e). 
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On the other hand, recent relativistic numerical simula- 
tions of spinning black hole binaries (Herrmann et al. 2007, 
Schnittman k, Buonanno 2007) suggest that Vr increases lin- 
early with the black hole spin parameter a , where ^ a ^ 1, 
and in the case of highly spinning black holes (a > 0.8) 
the magnitude of the kick suffered by the remnant could be 
of the order of a few thousand km/s (Tichy & Marronetti 
2007). Campanelli et al. (2007) report values of Vr as high 
as ~ 4000 km/s for equal mass binaries, if both spins lie 
in the binary orbital plane. Such a kick is sufficient to eject 
the remnant not only from a dwarf galgixy, where the escape 
velocity is ~ 300 km/s, but even from the center of a giant 
elliptical, for which the escape velocity can reach 2000 km/s. 

Though it is likely that MBHs acquire high spins (e.g. 
Volonteri et al. 2005) during their accretion history, the im- 
pact of the resulting recoil on the MBH cissembly has never 
been studied in details so far. If extreme recoil is indeed 
the rule, the ejection of a large fraction of MBHs formed 
through the coalescence of a binary systems can cause a sig- 
nificant drop in the number of expected coalescing events on 
the way of MBH assembly. Volonteri (2007) recently showed 
that current assembly models are able to reproduce the ma- 
jor observational constraints even if the extreme recoil pre- 
scription by Campanelli et al. (2007) is taken into account. 
However, high kick velocities could seriously affect the ex- 
pected number counts predicted for LISA, since the ejec- 
tion of remnants by their host halos would avoid subsequent 
MBHB formation. 

In this letter we estimate a robust lower limit for the 
predicted number of LISA sources. We use the Montecarlo 
realizations of the merger history performed by Volonteri 
(2007) to show that even in the worse (for GW observa- 
tions) case scenario in which during each merger the two 
MBH spins are counter-aligned in the MBHB orbital plane 
and extreme recoil is at work, current MBH assembly mod- 
els predict that at least ten sources will be dotoctablc by 
LISA. In practice, the lower limit of the expected rmmber 
of LISA sources does not substantially drop with respect 
to models employing non-spinning MBH recoil prescriptions 
(e.g. Volonteri Haardt & Madau 2003). 



Volonteri (2007) that are representative of these two classes 
of MBH assembly scenarios: the VHM and the BVRlf mod- 
els. In the VHM model, representative of the light seed sce- 
narios, (Sesana et al. 2007 for details) seed MBHs form with 
masses rrisccd ~ fewxlO^ M0, in halos collapsing at 2; = 20 
from rare 3.5-0" peaks of the primordial density field. In the 
BVRlf model, representative of the heavy seed scenarios, 
(Sesana et al. 2007 for details), black hole seeds form in 
halos subject to runaway gravitational instabilities, via the 
so-called "bars within bars" mechanism (Shlosniaii, Frank 
& Begelman 1989). MBH seed formation is assumed to be 
efficient only in metal free halos with virial temperatures 
Tvir^W^K, leading to a population of massive seed black 
holes with mseed ~ fewxlO^ Mq. 

The subsequent MBH evolution relies only on a few 
simple assumptions. Nuclear activity is triggered by halo 
mergers: in each major merger the more massive hole ac- 
cretes gas until its mass scales with the fifth power of the 
circular velocity of the host halo, normalized to reproduce 
the observed local correlation between MBH mass and ve- 
locity dispersion (mBH — cr, relation). MBHB coalescence is 
assumed to occur efficiently following halo mergers. 

For both the VHM and the BVRlf models, we consider 
two cases that bound the effect of recoil in the assembly of 
MBHs and, as a consequence, LISA events: (i) no gravita- 
tional recoil takes place and (ii) maximal gravitational recoil 
is associated to every MBHB merger, using the model by 
Volonteri (2007), which is based on the estimates reported 
by Campanelli at al. (2007) . For the latter we use the merger 
tree realizations presented in Volonteri (2007). The model 
takes into account consistently for the cosmic evolution of 
the mass ratio distribution of merging binaries and of their 
spin parameters (see discussion in Volonteri 2007). In each 
single merger, the mass ratio and the MBH spin magni- 
tudes are therefore fixed by the merger hierarchy; the spin 
orientations are instead chosen so as to maximize the recoil. 
MBH spins are assumed to initially lie in the binary orbital 
plane, counter-aligned one to each other. The recoil velocity 
is then evaluated according to equation 1 of Campanelli et 
al. (2007), that in this case simplifies as: 



2 MODELS OF BLACK HOLE FORMATION 

In the hierarchical assembly framework, MBHs form grow- 
ing trough mergers and accretion from seed black holes at 
high redshift. There are two main scenarios for MBH assem- 
bly, namely the light seed and the heavy seed models. In the 
light seed models, seed MBHs typically form with masses 
JTisocd ~ fewxlO^Me, in halos collapsing at 2 ~ 20, and 
are thought to be the end-product of the first generation 
of stars (Madau & Rees 2001). In the heavy seed models, 
black hole seeds form already massive (10* - 10^ M©) fi-om 
the low angular momentum tail of gas in protogalaxies at 
high redshifts. The angular momentum distribution of the 
gas in early-forming halos can be determined by means of 
cosmological A/'-body simulations (Bullock et al. 2001); halos 
with low spin parameters are prone to global dynamical in- 
stabilities, loading to the formation of a massive seed black 
hole (Koushiappas et al. 2004, Begelman Volonteri & Rees 
2006, Lodato & Natarajan 2006). 

We focus here on the two specific models discussed in 
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Here ^ = 1.2 x 10* km/s, B = -0.93, if = 6 x 10* km/s, oi 
and 02 are the magnitudes of the spin parameters of the two 

holes, ey is a unit vector in the binary orbital plane and e± 
defines the direction perpendicular to the orbital plane. The 
component of Vr along the e± direction depends sinusoidally 

upon the angle between the MBH spins and their initial 
linear momenta. To get the maximum recoil we set O = Oo 
(~ 0.184). 

We would like to emphasize that the prescription that 
we have chosen for (ii), and whose main features we have just 
summarized is the least favourable for gravitational wave ob- 
servations and (probably) unlikely to occur in these extreme 
circumstances during MBH assembly (Bogdanovic Reynolds 
& Miller 2007). 



Gravitational waves from massive black hole binaries 3 



3 GRAVITATIONAL WAVE SIGNAL 



Full discussion of the GW signal produced by an inspi- 
raling MBHB can be found in Sesana et al. 2005, along 
with all the relevant references. Here we just recall that a 
MBH binary at (comoving) distance r{z) with chirp mass 
M = m'l^^m2^^ / {mi + 01,2)^''^ generates a GW signal with 
a characteristic strain given by (Sesana et al. 2005): 
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An inspiraling binary is then detected if the signal-to-noise 
ratio (S/N) integrated over the observation is larger than a 
given detection threshold, where the optimal S/N is given 
by (Flanagan & Hughes 1998) 
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Here, / — fr/{l + z) is the (observed) frequency emitted at 
time f = of the observation, and the integral is performed 
over the frequency interval spanned by the shifting binary 
during the observational time. Finally, ftrms = •y/5/S'h(/) is 
the effective rms noise of the instrument; Sh{f) is the one- 
sided noise spectral density, and the factor \/5 takes into ac- 
count for the random directions and orientation of the wave; 
hims is obtained by adding in the instrumental noise contri- 
bution (given by e.g. the Larson's online sensitivity curve 
generator http: / /www.srl. caltech.edu / ~shane / sensitivity ) , 
and the confusion noise from unresolved galactic (Nelemans 
et al. 2001) and extragalactic (Farmer & Phinney 2003) 
WD- WD binaries. Notice that extreme mass-ratio inspirals 
(EMRI) could also contribute to the confusion noise in the 
mHz frequency range (Barack & Cutler 2004). 



4 RESULTS 

4.1 Coalescence rates 

Figure [T] shows the number of MBH binary coalescences per 
unit logA^ per unit observed year, dN/dlogMdt, predicted 
by the two models that we have considered, for both cases 
where recoil is neglected and extreme recoil is taken into 
account. Each panel shows the rates for different redshift 
intervals. Note that when extreme recoil in included, the 
rate predicted by the BVRlf model at any redshift is only 
marginally affected, while the VHM model is more sensi- 
tive to the GW recoil: at 2; > 15, GW kicks do not affect 
the coalescence rate; on the contrary, eX z < 15, the rate 
drops by a factor of ~ 3 for A^>10^ M©, if extreme kicks 
are included in the evolution. This is related to the fraction 
of seeds that experience multiple coalescences during the 
MBH assembly history. We can schematically think of the 
assembly history as a sequence of coalescence rounds, as also 
recently suggested by Schnittman (2007). After each round 
extreme recoil depletes a large fractions of remnants, and 
the relative importance of each subsequent round drops ac- 
cordingly. In the VHM model, about 65% of the remnants of 
the first round will undergo a second round of coalescences, 
so the second round has an important relative weight in the 
computation of the total rate. When extreme recoil is taken 
into account, a large fraction of the first round remnants is 
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Figure 1. Number of MBHB coalescences per observed year at 
2 = 0, per unit log chirp mass, in different redshift intervals. 
Solid lines: GW recoil neglected; dashed lines: extreme GW recoil 
included. Thick lines: VHM model; thin lines: BVRlf model. 



ejected from their hosting halos. We find that the effective 
fraction of remnants that can experience a second coales- 
cence drops to ~ 30%. This is the reason why the number of 
coalescences involving light black holes (A^ < 10'' Mq) does 
not drop at any redshift, while the number of coalescences 
involving more massive binaries drops by a factor ~ 3. In 
the BVRlf scenario seeds are rarer, and the fraction of first 
coalescence remnants that participate to the second round 
is around 25%; switching on the extreme recoil has a signif- 
icantly smaller impact on the global rate in this case. More- 
over, in this model seeds are more massive and the bulk of 
merging events happens at lower redshift, where the hosting 
halo potential wells are deeper and consequently larger kicks 
are needed to eject the coalescence remnants. As a matter 
of fact, the seed abundance sets the mean number of ma- 
jor mergers that a seed is expected to undergo during the 
cosmic history, and this basically sets the ability of extreme 
kicks to reduce the coalescence rate. 



4.2 LISA detection rate 

We now discuss how the number of GW sources detectable 
by LISA is influenced by extreme GW recoils. To facili- 
tate the comparison with our previous works, all the results 
shown here assume an observation time of 3 years, a sharp 
low-frequency wall at lO"** Hz in the instrumental sensi- 
tivity (see Sesana et al. 2007), and a detection threshold 
S/N = 5 (see equation |3]) ; the confusion noise includes only 
galactic and extragalactic white dwarfs and ignores a possi- 
ble contribution from EMRIs (Barack & Cutler, 2004). At 
the end of this section, we will briefly discuss the impact of 
the former assumptions on the number of detectable sources. 
Figure [2] shows the redshift distribution of MBHBs detected 
by LISA. The effect of extreme GW recoils on the source 
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Figure 2. Redshift distribution of MBHBs resolved with S/N > 
5 by LISA in a 3-year mission. Line style as in figure [T] The top- 
right corner label lists the total number of expected detections. 
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Figure 3. Chirp mass function of MBHBs resolved with S/N > 5 
by LISA in a 3-year mission. Line style as in figure [l] All curves 
are normalized such as the integral in dlog(jVf) gives the number 
of detected events. 



number counts drastically depends on the abundance and 
nature of the seeds, along the lines discussed in the previ- 
ous section. In the VHM model, the number of detectable 
sources drops by a factor ~ 60%, and the number of the 
potential LISA detections is reduced from « 140, if the re- 
coil is neglected, to « 60, if extreme recoil is included. Vice 
versa, the detection rate predicted by the BVRlf model is 



Figure 4. Impact of LISA sensitivity details on the number of 
detected sources. Solid lines: S/R = 5 and sensitivity cut-off at 
/ = 10"'* Hz; long-dashed lines: S/R = 8 and sensitivity cut-off 
at / = 10""* Hz; short-dashed lines: S/R = 8, sensitivity cut-off 
at / = IQ--* Hz and EMRI confusion noise added; dotted lines: 
S/R = 5 and sensitivity cut-off at / = 10~^ Hz. Thick lines are 
for VHM model, thin lines are for the BVRlf model. The number 
of detected MBHBs in a 3 year observation, under the different 
detection assumptions, are also listed. Extreme recoil is assumed. 



only weakly affected by the extreme recoil prescription, and 
it drops by about 15% (from 40 to 34 events in 3 years of 
observation). Note that though the overall number of coa- 
lescences in the VHM model decreases only by about 25% 
when extreme recoil is considered, the number of LISA de- 
tections is reduced by a much larger factor. This is because 
if the seeds are light, LISA can not detect the bulk of the 
first coalescences of light binaries happening at high red- 
shift, that are responsible for the major contribution to the 
coalescence rate and are not affected by the recoil. LISA can 
observe later events, involving more massive binaries, that 
are largely suppressed by the MBH depopulation due to ex- 
treme GW kicks. In the BVRlf model, on the other hand, 
seeds are more massive, and the second coalescence round is 
less important; in this case, the LISA sensitivity is sufficient 
to observe almost all the first coalescences, and the number 
of detections is only mildly reduced. As the kicks affect the 
merger rate starting from the second round, its signature 
consists in a slight decrease of the mean chirp mass of the 
detected binaries, see figure |31 

We emphasize here two aspects (i) at time it is not clear 
if LISA will be able to shed light on the importance of re- 
coil in MBH assembly, even in this extreme case, since the 
uncertainty introduced in the number counts is at most of 
a factor of ~ 3, comparable with uncertainties due to our 
ignorance in the MBH accretion history and in the detailed 
dynamics of MBHBs (see, e.g., discussion in Sesana et al. 
2007); (ii) on the other hand, this fact confirm that MBHBs 
are LISA safe targets; since extreme recoil effects increase 
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with the seed abundance, we expect the drop in the detec- 
tions to be more significant for those scenarios that predict 
a larger number of sources. In figure |4] we show how differ- 
ent assumptions on the detection threshold, the instrumen- 
tal noise below 10~* and the confusion noise from EMRI 
affect the LISA detection rates. If seeds are massive, the re- 
sults shown in figure [2] are hardly affected. If seeds are light, 
EMRI confusion noise and a more conservative detection 
threshold, say S/N — 8, can halve the number of sources 
detected by LISA. For both scenarios, extending the LISA 
sensitivity window below 10~* has also minimal effect on 
the number of detections 



5 DISCUSSION 

Here we have considered two specific MBH assembly models, 
representative of two different MBH seed formation scenar- 
ios. However our findings can be considered, at least quali- 
tatively, valid in general. Given the size and the abundance 
of the seeds, our 'coalescence round' picture depends on the 
details of the models. For example in the VHM model we 
checked that by changing the accretion prescription (see 
Volonteri, Salvaterra & Haardt 2006) the total number of 
events would change by a factor of two (note that the ac- 
cretion prescription considered in the models described in 
the previous section gives the minimum number of coales- 
cences); however the relative weights of the different coales- 
cence rounds do not change significantly. So we can safely 
conclude that a decrease ^50% in the expected LISA sources 
should be a general trend for all those models in which 
the MBH assembly starts from light seeds at high redshift. 
In this class of models the number of predicted coalescing 
events is so high (>100 yr~^) that at least a few tens of 
MBHBs should be guaranteed LISA sources. On the other 
hand, extreme recoils should not be an issue at all for LISA 
if the MBH seeds are massive and/or rare. We remark here 
that in the BVRlf model we assumed MBH seed formation 
to be efficient only in metal free halos with virial temper- 
atures rvir>IO''K, i.e, we have considered atomic hydrogen 
to be the only coolant. Assuming efficient molecular hydro- 
gen gas cooling (e.g. Koushiappas, Bullock & Dekel 2004) 
the number of seed MBHs increases by an order of magni- 
tude (and being the seeds massive, LISA would be able to 
detect the first coalescence round), and the GW kick would 
not be an issue at all. Relying on this results, the estimate 
of ~ 10 detections in three years predicted by the BVRhf 
model described in Sesana et al. 2007 does not change un- 
der the assumption of extreme recoils (seeds are heavy and 
rare), and can be considered a robust LISA detection lower 
limit. To conclude, in Sesana et al. 2007 we explored differ- 
ent MBH assembly scenarios to quantify the imprint of the 
MBH seed prescription on the LISA data stream. Motivated 
by recent studies on extreme GW recoils, we have quantified 
in this letter their impact on the MBHB coalescence and on 
the LISA detection rate, confirming that the detection of at 
least ~ 10 coalescing binaries in a 3 year mission is a robust 
prediction even considering extreme GW recoils. 

I would like to thank M. Volonteri for having provided 
the outputs of the merger trees used in this study and A. 
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on the manuscript. 
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